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Abstract 

We study the evolution of massive scalar field in the spacetime geometry of Kehagias-Sfetsos(KS) 
£>> black hole in deformed Hofava-Lifshitz(HL) gravity by numerical analysis. We find that the signa- 

ture of HL theory is encoded in the quasinormal mode(QNM) phase of the evolution of field. The 
QNM phase in the evolution process lasts for a longer time in HL theory. QNMs involved in the 

CM 

evolution of massive field are calculated and find that they have a higher oscillation frequency and 

CJ . 

a lower damping rate than the Schwarzschild spacetime case. We also study the relaxation of field 
in the intermediate and asymptotic range and verified that behaviors of field in these phases are 



independent of the HL parameter and is identical to the Schwarzschild case. 



Electronic address: nijovarghese@cusat.ac.in, vck@cusat.ac.in 



I. INTRODUCTION 



One of the major open challenges facing theoretical physics is to accommodate general 
relativity (GR) in the framework of quantum field theory The main problem one encounters 
is the non-renormalizablity of gravity at the perturbative level. Horava[lJ suggested a power- 
counting renormalizable theory of gravity in 3+1 dimensions inspired from the Lifshitz model 
now dubbed as Horava-Lifshitz gravity. It assumes a Lifshitz-like anisotropic scaling between 
space and time at short distances, characterized by a dynamical critical exponent z = 3 and 
thus breaking the Lorentz invariance. While in the IR limit it flows to z = 1, retrieving 
the Einstein's GR. Even though many of the fundamental questions have not yet been 
answered the seemingly innovative ideas of HL theory ignited an interest in the topic 
and a number of works has been reported based on it. Since the theory has the same 
Newtonian and post Newtonian corrections as those of GR, systems of strong gravity, like 
black holes, are needed to get observable deviation from the standard GR. Various black 
hole solutions are found in HL theory of which one with asymptotically flat Minkowski 



spacetime is the KS black hole given in 



24| by applying de 



ormation in the original theory. 
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37|. 



Various aspects of KS black hole were explored in the past 

The study of response of the black hole to external perturbations is the best way to get a 
clear picture of the properties of these exotic objects. As it is well established, the evolution 
of perturbations in black hole spacetimes involves three stages 461] . The first one is an initial 
response containing the information of the particular form of the original wave field followed 
by a region dominated by damped oscillation of the field called quasinormal modes, which 
depends entirely on the background black hole spacetimes. The QNMs play an important 
role in almost all astrophysical processes involving black holes and the spectra of QNMs may 
be different in various theories of gravity and would help us to distinguish these theories. 
The late time tail stage dominates in the final phase of evolution. 

Evolution of different fields and the associated QNMs were studied in HL black holes in 



the frequency domain 



38 



44l | . All these studies on the field evolution 



and in time domain 

around black holes in HL theory has so far been restricted to the massless fields. Evolution 
of massive field is a more intriguing subject of study. The spectrum of entropy/area is 



discussed in the viewpoint of QNMs of massive scalar field in 



45| . It is well understood in 



standard GR that massive field behaves differently from the massless one. It was confirmed 
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that the massive fields decay more slowly than massless one 49J. The late-time behavior of 
the field that follows the QNM stage of evolution is also an interesting topic of study since 
it reveals the actual physical mechanism by which a perturbed black hole sheds its hairs. It 



was first demonstrated by Price 61( that in the background of the Schwarzschild spacetime 
the massless neutral field at late time dies off as \I/ ~ £~( 2 ^+ 2 ) or ^ ~ £-( 2 ^+3) depending 
on the initial conditions and the multipole order t. In contrast with the massless fields, 
massive fields have an oscillatory inverse power-law behavior, with \& ~ t~^ +3 / 2 ),sm(mt) at 
intermediate late times Mm <C mt <C l/(mM) 2 [48]. But it was shown analytically that 
in the asymptotic late times (mt ^> l/(mM) 2 ) another pattern of oscillatory tail of the 
form ^ ~ t~( 5//6 ) sin(mt) dominates 50| and it was numerically verified for various space 



times and fields 



52H55j . All these studies argued that late time relaxation does not have any 



relation to the space time parameters. Inspired from the hints that QNMs show dependence 
on spacetime parameter in three dimensional AdS|56| and dS 57] space times, a detailed 
numerical study of relaxation process in RN space time was done in 58]. They showed 
that for Mm <C 1 relaxation depends only on the field parameters, but when Mm ^> 1 
spacetime parameters affect the relaxation and found that for a Schwarzschild black hole 
bigger the black hole mass is, the faster the perturbation decays. The purpose of this paper 
is to investigate the different stages of evolution of a massive scalar field around black hole 
in HL theory. 

The paper is organized as follows. In SecJTT] we derive the wave equation for massive 
scalar field around KS black hole and the numerical method employed to study the time 
evolution is explained and in Sec JIIII the evolution of a massless field is studied using the 
method of time domain integration. The massive field perturbation is studied in Sec JIVI and 
the results are summarized in SecJVl 



II. SCALAR FIELD AROUND KS BLACK HOLE 



In |24j Kehagias and Sfetsos obtained a spherically symmetric black hole solution by 
deforming HL gravity that allows asymptotically flat spacetimes with the addition of a term 
proportional to the Ricci scalar of three-geometry, fj^R^ 3 ', with the cosmological constant 
Aiy — > 0. This will not alter the UV properties of the theory but it does the IR ones 
leading to Minkowski vacuum analogous to Schwarzschild spacetime in GR. The external 
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gravitational field of such a black hole has the form 24j , 

ds 2 = -N(r) 2 dt 2 + f(r)- l dr 2 + rW, 
where f(r) has the form, 

N(r) 2 = /(r) 
The event horizons are at, 



2(r 2 - 2Mr + a) 



- 2 + 2a + Vr A + 8aMr' 



(1) 



(2) 



r± = M± VM 2 - a. (3) 

When the HL parameter a = the solution reduces to the Schwarzschild spacetime. We 
will now consider the evolution of scalar field in this spacetime. Scalar field evolution is 
governed by the Klein-Gordon equation, 



'-9 



^(v^^)$-m 2 $ = 0, 



(4) 
and 



Resolving the scalar field into spherical harmonics, <&(£, r, 9, 4>) = -^i(t,r)Yi m (9, 
employing the tortoise coordinate defined by dr* = jdr, the above equation can be reduced 
to the form, 



where the effective potential Vi(r) is given by 



■V e {r)V e {t,r) = Q, 



(5) 



VAr) 



f(r) 



+ 1) ldf(r) 



r dr 



m 



(6) 



The behavior of the effective potential Vg{r) is plotted in FigJTJ In FigJIta) the dependence 
of the potential with a is shown for massless and massive(m = 0.3) fields. For a given mass, 
as the parameter a increases the peak of the potential increases where as the asymptotic 
region of potential is unaffected. It is clear from Figljjb) that as the mass of the field 
increases the hight of the potential increases and its asymptotic value raises as m 2 . The 
barrier nature of potential will be spoiled for high field masses. 

The complex nature of the potential makes it difficult to obtain the exact solutions of 
Eq([5]) and we have to tackle the problem by numerical methods. The wave equation Eq(j5]) 
can be recasted in the null coordinates u = t — r* and v = t + r* as, 

. d 2 



dudv 



ty(u,v) = V(u,v)^(u,v). 



(7) 



|* a -j 2 4 6 8 10 12 14 ^) 5 10 15 20 25 30 

FIG. 1: Profiles of effective potential Ve(r) for 1 = 2. (a) Vi(r) for a = 0.2(solid) and a = 
0.8(dotted). (b) Vp(r) with different field masses for a = 0.4. 

The above equation is invariant under the resealing 

r — 7- ar, t — )■ at, M — > aM, u — >■ u/a, a — > a 2 a (8) 
for some positive constant a, such that Mm is invariant. We numerically integrate this 



equation using the difference scheme 



■^ N = ^ W + ^ E -^ S - ^ V {S){^ W + * E ) + 0(h A ) (9) 

o 

The integration is performed on an uniformly spaced grid with points N(u + h,v + h), 
W(u + h,v), E(u, v + h) and S(u, v) forming a null rectangle with an overall grid scale factor 
of h. Since the late-time behavior of the wave function is found to be insensitive to the 
initial data, we set ip(u,v = 0) = and a Gaussian profile ty(u = 0, v ) = Aexp — ^ • 
In all our calculations we set the initial Gaussian with width a = 3 centered at v$ = 10. Due 
to the linearity of equation one has the freedom to choose the amplitude of the initial wave 
A = 1. To proceed the integration on the above numerical scheme one has to find the value 
of the potential at r(r*) = r((v — u)/2) at each step. For this, we numerically integrate the 



equation for the tortoise coordinate using the Runge-Kutta method 
interpolation, obtained r(r*) at each step. 



601 ] and by cubic spline 
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III. EVOLUTION OF MASSLESS FIELD 



Before studying the massive case, we first examine the massless scalar field(m = 0) 
evolution around KS black hole. In FigfSJ^a) time evolution of a typical 1 = 2 field in HL 
theory is compared with the corresponding Schwarzschild case. We can see that the QNM 
phase extends for a longer time in HL theory. The late-time tail starts at a later time t ~ 415 
for KS black hole with a = 0.8 whereas it is at t ~ 320 for the Schwarzschild case. It is 
also clear from the figure that the oscillation frequency and the damping time have a higher 
values in HL theory. 

After the QNM stage the massless field dies off as an inverse power of time \I/ ~ t~( 2£+3 \ 
FigfSJJb) shows the late-time behavior of wave function for different values of a, with multi- 
pole index i = 2. It is clear from the figure that the late-time behavior massless field decay 
is independent of a and is same as for the Schwarzschild black hole. 




2.5 3 3.5 4 4.5 5 5.5 6 6.5 (J>) 51 5 - 8 5 - 9 6 6 - 1 62 6 - 3 64 65 



FIG. 2: Temporal evolution of massless scalar field in log-log scale with 1 = 2. (a)Wave function 
in KS spacetime with a = 0.8 in comparison with that in Schwarzschild spacetime(o = 0). (6)Late 
time behavior of the field for deferent values of a. 

In FigJHJ field evolution for different multipole indices are explored for a = 0.4. The 
field falls off as tp ~ t~ 3 - 14 , t~ 5 - 13 , t~ 7A2 and t~ 913 for t = 0, 1,2 and 3 respectively with the 
predicted powers —3, —5, —7 and —9 respectively. 
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2.5 3 3.5 4 4.5 5 5.5 6 6.5 
t 



2.5 3 3.5 4 4.5 5 5.5 6 6.5 
t 



FIG. 3: Decay of field having different £. The late-time tail dies off with power-law exponents 
-3.14, -5.13, -7.12 and-9.13 for i = 0, 1,2 and 4 respectively. 

IV. EVOLUTION OF MASSIVE FIELD 

In this section we consider the massive field evolution in the spacetime of KS black hole. 
We will compare the results with the existing results in standard GR. Here in the first series 
of experiments we set the mass of the black hole, M = 1 and choose the the field mass m such 
that Mm < 1. In FigJ3|a)the wave function for massive scalar field(m = 0.1) is plotted in 
comparison with the corresponding Schwarzschild case for i = 2. After the prompt response 
in the beginning the quasinormal ringing starts. We can see a clear difference between the 
two cases in the QNM region that the QNM phase ends at a later time in HL theory and 
has a lower damping rate. FigJD^b) shows the variation of QNM with the field mass m for 
a = 0.4. As the mass of the field increases the QNM phase shrinks in time. 

To see the effects of the HL parameter a and field mass m in QNM phase, we calculate 
the exact values of QNMs from the numerically integrated data in the time domain by a 



nonlinear x 2 fitting. We also use the third order WKB method 63[, an efficient semianalytic 
method originally developed by Schutz and Willjo^l for the lowest order. The method is 
found to be accurate for low lying modes and can be used to explore the QNM behavior of 
black holes for field with very low mass. The calculated values are given in Table |T] and Table 
HT1 QNMs in HL theory have a higher oscillation frequency Re(cu) and a lower damping rate 
\Im(u)\ than the Schwarzschild spacetime case. Also the Re(co) and \Im(u)\ found to be 
decreasing with the increase of the HL parameter a. WKB method gives verse results for 
higher field masses since the barrier nature of potential gets spoiled at these mass ranges. 
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2.5 3 3.5 4 4.5 5 5.5 6 6.5 
<»> ln(t) 



4 5 
ln(t) 



FIG. 4: Time evolution of massive scalar field for I = 2 mode. (a)Field of mass, m = 0.1 in KS 
spacetime with a = 0.8(top curve) in comparison with the corresponding case in Schwarzschild 
spacetime(bottom curve). (fe)Field with different masses for a = 0.4, curves from bottom to top is 
for m = 0.01, 0.05, 0.1 and 0.2. 

So there should be a discrepancy between the QNMs evaluated by time domain and WKB 
methods at high field mass. 







WKB 


Time domain 


a 


Re{uj) 


Im(uj) 


Re(oj) 


Im{ui) 





0.48637 


-0.09572 


0.48552 


-0.09573 


0.2 


0.49388 


-0.09253 


0.49377 


-0.09256 


0.4 


0.50231 


-0.08884 


0.50159 


-0.08888 


0.6 


0.51188 


-0.08438 


0.51142 


-0.08437 


0.8 


0.52293 


-0.07868 


0.52143 


-0.07893 


1 


0.53587 


-0.07069 


0.53569 


-0.07106 



TABLE I: Fundamental(n = 0) QNM frequencies of massive(m = 0.1) scalar field for i = 2, 
calculated using WKB and numerical integration data 

The QNM phase is followed by the phase of late time tail behavior of the field decay. First 
we consider the late time tails in intermediate range, Mm mt (mM) 2 . The decay of 
field along black hole outer horizon H + (approximated by the null surface u max = 0.5 x 10 4 ) 
and future time like infinity i + (approximated on fixed radius r* = 50M) are evaluated with 
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WKB 


Time domain 


Tfl 


Re(oj) 




1 LC- \ UJ j 


/ 7T1 1 t • ) 1 
J / / fc 1 UJ J 


n 


D 49942 


-0 08969 


u.^yoou 


-u.uoyu i 


0.02 


49953 


-0 08966 


D /LQQ7Q 


-u.uoyuy 


0.05 


50014 


-0 08948 




-u.uoy^o 


0.07 


50083 


-0 08928 


U.OUXtJU 




0.1 


50231 


-0 08883 


0.50159 


-0.08875 


0.2 


0.51101 


-0.08622 


0.51083 


-0.08701 


0.3 


0.52571 


-0.08169 


0.52737 


-0.08029 


0.4 


0.54665 


-0.07501 


0.53933 


-0.06758 


0.5 


0.57423 


-0.06568 


0.54165 


-0.05451 



TABLE II: QNM frequencies afor different mass of the scalar field m. QNMs calculated using 
WKB and numerical integration data (£ = 2 and a = 0.4) 

the initial field profile parameters vq = 50 and cr 2 = 2. The result for a = 0.4 with m = 0.01 
and £ = is shown in Figj^a). After the QNM phase the amplitude of the oscillatory field 
decays as a power law with exponent -1.49 along I + and H + with a period of oscillation 
T = 314.8. We have analyzed for different values of HL parameter a and found that the 
tail behavior is almost independent of a and the decay is identical to the Schwarzschild case 
according to the form, 

^ _ r («+3/2) sin ( mt ) ( 10 ) 

This can be expected since the late-time tail is originated by the backscattering by the 
effective potential in the asymptotic region and as it is clear from Fig^TJ that the parameter 
a has no effect on this asymptotic region of the potential but can only change its shape near 
the peek. 

The field decay for different multipole indices is plotted in Figj5|b) with m = 0.01. The 
decaying tail is found to have a period of oscillations T = 314.7 ± 0.03 and the power-law 
exponents —1.49,-2.5,-3.51 and —4.51 for I = 0,1,2, and 3 respectively. An excellent 
agrement with the decay rate of the form £~^+ 3 / 2 ) can be seen. We find that the frequency 



9 




FIG. 5: Decay of massive field with m = 0.01. (a) Field along outer horizon H + (top curve) 
and future time like infinity i+ (bottom curve). The amplitude of field decay as power law with 
exponent —1.49 and period of oscillation T = 314.8(6) Field along i + for different t. The field dies 
off as power-law with exponents —1.49,-2.5,-3.51 and —4.51 for t = 0,1,2, and 3 respectively. 
The period of oscillation T = 314.7 ± 0.03. 

and damping rate of the oscillatory tails are independent of the parameter a and follow the 
Schwarzschild case. Thus Eqf TTOj) will be the form of the intermediate late-time behavior of 
massive field decay in HL theory also. 




FIG. 6: Late time behavior of field for a = 0.4. (a)Maxima of oscillation for I = 0, from bottom 
to top m = 0.02,0.05,0.07. Solid line has a slop —1.5. (fe)Asymptotic regime of field decay with 
m = 1 for I = 0(top) and I = 2(bottom). Straight lines have a slop —5/6. 

But in the asymptotic late time the decay of field does not follow pattern given by Eq ffTUl) . 
Another pattern of oscillatory tail dominates in the asymptotic late times mt ^> l/(Mm) 2 . 
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This can be seen from Figfj^a), where the maxima of oscillations are shown for a = 0.4, 1 = 
and m = 0.02,0.05,0.07. The deviation from straight line with slop —1.5 is visible in the 
asymptotic late time. When mt ^> l/(Mm) 2 smaller the Mm the later the asymptotic tail 
starts. In FigJ^b) asymptotic region is shown for m — 1 and I = and 2. We find that 
the asymptotic tail is independent of the multipole order and follows the \1/ ~ t~ 5 ^ 6 sin(mt) 
form as in the standard GR. We have analyzed the asymptotic tail for different values of a 
and find that this regime is also independent of the HL parameter a. 




ln(l) 



FIG. 7: Decay of field with different black hole mass for m = 0.01, a = 0.4 and I = 0. For 
M = 0.5,1 and 1.5(curves from top to bottom) the decay rates of field are —1.458,-1.459 and 
—1.462 respectively. 



Now we study the role of the black hole parameter on the decay of field for Mm <C 1. 
The numerical results for different black hole mass with field mass m = 0.01 are shown in 
FigJTl The field falls off with power law exponents —1.458, —1.459 and M = 0.5, 1 and 1.5 
respectively and we can conclude that for Mm 1 the decay rate is independent of black 
hole mass. We are not going for the Mm ^> 1 case of computationally expensive part of 
study. For Mm ^> 1 range we have good reasons to believe that the relaxation will show 
similar results of the Schwarzschild case as shown in 58(, since we have already shown that 



the parameter a has no effect on the asymptotic region of the potential and the late time 
relaxation is independent of the HL parameter. 
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V. SUMMARY 



We have studied the time evolution of massive scalar field in the spacetime geometry of 
KS black hole in deformed HL gravity by time domain integration. We find a noticeable 
deviation of the evolution behavior of massive field in the ringdown region from the standard 
Schwarzschild black hole. QNMs involved in the evolution of massive field are calculated 
from the time domain integration data and third order WKB method. Massive QNMs in HL 
theory have a higher oscillation frequency and a lower damping rate than the Schwarzschild 
spacetime case. The ringdown frequency and damping rate also decrease with the increase 
of the HL parameter a. As the field mass increases the QNM phase squeezes to a smaller 
time interval. However the late time evolution of massive field fails to show any distinction 
from the Schwarzschild case. In the intermediate range the field decays as t~^+ 3 / 2 ) sin(mt), 
but in the asymptotic late time the decay is dominated by t~ 5 / 6 sin(mt) tail. In the range 
mt ^> l/(Mm) 2 , smaller the value Mm the later the asymptotic tail starts. We also shown 
that for Mm <C 1 the relaxation is independent of space time parameters. Since the HL 
parameter a has no effect on the asymptotic region of the potential the late time relaxation 
is independent of the HL parameter. 
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